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Abstract 

We study the maximal entropy per unit generator of push-point mapping classes 
on the punctured disk. Our work is motivated by fluid mixing by rods in a pla- 
nar domain. If a single rod moves among A^-fixed obstacles, the resulting fluid 
diffeomorphism is in the push-point mapping class associated with the loop in 
7ri{D^ — {N points}) traversed by the single stirrer. The collection of motions in each 
of which the stirrer goes around a single obstacle generate the group of push-point 
mapping classes, and the entropy efliciency with respect to these generators gives a 
topological measure of the mixing per unit energy expenditure of the mapping class. 
We give lower and upper bounds for Eff(A^), the maximal efliciency in the presence 
of obstacles, and prove that Efl(A^) log(3) as ^ oc. For the lower bound 
we compute the entropy efliciency of a speciflc push-point protocol, HSPat, which 
we conjecture achieves the maximum. The entropy computation uses the action on 
chains in a Z-covering space of the punctured disk which is designed for push-point 
protocols. For the upper bound we estimate the exponential growth rate of the ac- 
tion of the push-point mapping classes on the fundamental group of the punctured 
disk using a collection of incidence matrices and then computing the generalized 
spectral radius of the collection. 

1 Introduction 

The topological entropy h{f) of a map / is a standard measure of its dynamical complexity. 
The entropy of an isotopy class is the infinimum of all the entropies of the homeomorphisms 
in the class. If a surface isotopy class has positive entropy, a theorem of Thurston guarantees 
the existence of an (almost) canonical map in the isotopy class which realizes the infinimum 
and whose dynamics are described by a finite collection of mixing subshifts of finite type 
([Thu88], [FLP79]). Further, Handel's theorem shows that these dynamics are present up to 
semiconjugacy in every homeomorphism in the mapping class ([Han85], [Boy99]). 

The prototypical positive entropy mapping class is of pseudo- Anosov type and the Thurston 
representative in the class is a pseudo- Anosov map 0. In this case, A = exi[){h{(j))) is the 
dilation of the pseudo- Anosov map, an algebraic integer with a wide variety of geometric 
and topological interpretations. As a consequence of these interpretations it is natural to 
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seek the pseudo-Anosov map with least dilation within semi-groups of isotopy classes (see, 
for example, [HK06], [FLM09], or [LTll]). From a dynamics perspective, maximizing com- 
plexity is also of interest. Since entropy is multiplicative under iteration, h{f^) = nh{f)^ the 
entropy over a semigroup is always unbounded, and thus one needs a normalization. In many 
cases the semi-group of interest has a natural set of generators, for example, a collection of 
Dehn twists or motions of stirring rods as discussed below. In these cases following Thif- 
feault/Finn and Moussafir ([TF06], [Mou06], [FTll]) we consider the entropy efficiency of 
a mapping class which is the topological entropy per unit generator. The entropy efficiency 
of the semigroup is then the supremum of the efficiencies of its members. 

As a consequence of Thurston's theorem, the entropy of an isotopy class is equal to the 
exponential growth rate of words under iteration by the induced map on the fundamental 
group of the surface. One may also define this growth rate and thus the entropy of any 
automorphism of a finitely generated group G. Given a finite collection of automorphisms 
A — {ai, . . . , an} of the group, we define the entropy efficiency of the semigroup generated 
by A as the maximal entropy per unit generator (see §2.2). If G is a free Abelian group, 
the automorphisms are linear isomorphisms and the entropy efficiency of A is the log of the 
generalized spectral radius of the collection of linear transformations (Remark 2.1). The 
generalized and joint spectral radii are very useful and much-studied linear objects, and the 
entropy efficiency is their non-linear analog. 

Our main motivation for studying the entropy efficiency is fluid mixing via the motion 
of stirring rods in a planar domain ([BASOO], [FTll], [TF06], [Boy05]). In the simplest 
situation the mixing is accomplished by n-rods which after time-T return to their original 
positions. Such a stirring protocol is naturally described as a braid with n-strings. When 
the rods move through the fluid there is a resulting motion of the entire fluid region, and 
the evolution of a fluid particle from its initial position to its position after time-T deflnes 
the fluid motion diffeomorphism. The exact nature of this diffeomorphism depends on the 
physical properties of the fluid and the geometry of the fluid region and rod paths. However, 
as long as the fluid does not tear, all possible fluid motion diffeomorphisms arising from 
the same stirring protocol are in the same mapping class, and this class is the one usually 
associated with the braid (see §2.3). Thus if we can determine the entropy of this isotopy 
class we have a lower bound for the complexity of the evolution of any fluid motion resulting 
from the stirring protocol. The topological entropy has a particular physical signiflcance 
because it measures, among other things, the exponential growth rate of material lines in 
the fluid under the action of the protocol. Exponential stretching of material lines does not 
guarantee good mixing, but in many situations this rate of stretching is a good topological 
"flrst order" indicator. 

The other consideration in implementing a stirring protocol is its energy cost. Protocols 
which maximize the mixing per unit expenditure of energy could have great practical value, 
he energy required to implement a stirring protocol certainly depends on the geometry of 
the path of the stirrer, but since we are characterizing protocols just by their braids, we 
seek a topological measure of effort which just depends on the braid. Thus we pick some 
class of fundamental motions of the stirrers, for example, the switching of a pair of rods, and 
treat these as generators of a semigroup of mapping classes. The corresponding maximal 
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entropy efficiency gives a rough notion of the maximal mixing per unit expenditure of energy 
possible in the class of protocols. A protocol which realizes the maximum would be optimal 
in a precise mathematical sense, but perhaps not in any true physical sense, see [FTll]. 

The entropy efficiency of the full braid group with respect to the standard generators 
was studied by Finn and Thiffeault in [FTll]. They show that the maximal efficiency is 
bounded above by log[(l + \/5)/2]. Further, when n > 4, this bound is strict, and when 
n = 3 and 4 the bound is realized by (Jia2^ and erf ^0-20-3" ^a2, respectively. They also studied 
a generating set consisting of words of the form a^^^ . . . a^J^^ where for all j, \ij — > 1. 
For this class they obtained log(l + \/2) as an upper bound for the maximal efficiency. This 
class of generators models motions of the collection of rods in which the entire motion can 
be done simultaneously. 

Here we study the class of stirring protocols in which there are A^-fixed obstacles and a 
single stirrer moves. Thus the stirring protocols are determined by the path of the single 
stirrer in the A-punctured disk. This path can be considered an element of the fundamental 
group and so these stirring motions are called TTi-stirring protocols. Each vri-stirring proto- 
col determines a braid on (A' + 1) strings and an isotopy class on the (A' + l)-punctured 
disk. The mapping classes which result from vri-stirring protocols are called point-push in 
the mathematics literature ([Bir74], [KraSl], [DowlO], [FMar]). Thus the collection of tti- 
stirring protocols, denoted here as Vn^ is the subgroup of the braid group on (A + 1) strings 
corresponding to the point-push isotopy classes in the A^-punctured disk. To maintain the 
traditional mathematical terminology while incorporating the physical intuition and moti- 
vation of our problem, we will usually use point-push protocol to refer to elements of Vn- 

We use the analog of the standard generators of the fundamental group of the A^- 
punctured disk as our fundamental stirring motions. Thus for j = 1,...,A, let aj cor- 
respond to moving the stirrer counterclockwise around the j^^-obstacle (Figure 1), and we 
consider the entropy efficiency of a point-push protocol with respect to the generating set 
{ofi, . . . , ajsi} of Vn' The maximal entropy efficiency over Vn is denoted Eff(A'). Our main 
theorem is: 

Theorem 1.1 For N > 2 there are explicitly defined matrices H^^"^ and H^^"^ such that the 
maximal efficiency Eff(A') of a point-push protocol in the presence of N obstacles satisfies 



and so Eff (A) log(3) a5 A^ ^ 00. 

Note that the lower bound has no content when N — 2. The most efficient protocol with 
two obstacles is given in Theorem 4.7. 

' The proof of Theorem 1.1 follows directly from the lower bound given in Theorem 4.5 
and the upper bound given in Theorem 5.9. Both the upper and lower bound proceed by 
reducing the problem to the computation of the spectral radius of a specific product of 
generating matrices, but the reduction in the two cases is quite different. In both cases the 
required computations are tractable because of the simple forms of actions made possible by 
the adapted set of generators shown in Figure 2. 



log(3^ - 3A^ - 1) - log(A^) log(p(g(^))) 
N - N 



< Eff(A^) < 



log(p(^W)) ^ log(3^ - 2) 
A^ - A^ 
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Since the maximal efficiency is the supremum over the efficiencies of all point-push pro- 
tocols, the efficiency of any single protocol serves as a lower bound. In §4.1 we define a 
point-push protocol HSPat which numerically appears to give the maximal efficiency (Fig- 
ure 3). The invariant train track for HSPaa is rather simple, but the corresponding Markov 
transition matrix has very large entries. Thus the computation or even estimation of the 
dilation of HSPat as its spectral radius appears to be quite difficult. Instead we use homo- 
logical methods developed by Fried ([Fri86], [Fri83], [BB07]). In short, the idea is to lift the 
mapping class to finite covers where the dilation is computable as the spectral radius of the 
action on homology. This cover is found as the quotient of a Z-covering space whose homol- 
ogy is treated as a free Z[t^^]-module. The action of the mapping class is then a matrix over 
Z[t^^], and substitutions of t = exp{27rip/q) give pieces of the action of the mapping class 
on the homology of the corresponding Z^-covering space. 

As noted above, the entropy of a point-push protocol is equal to the exponential growth 
rate of words under iteration by the induced map on the fundamental group of the punctured 
disk, 7ri{DN^i). For each generator aj of Vn, we find the incidence matrix A^^^ of its action 
on 7Vi{Dn^i). The incidence matrix just counts the number of occurrence of a generator 
in the image of another generator. As also remarked above, the entropy efficiency is the 
nonlinear analog of the generalized spectral radius and in Proposition 5.1 we show that 
for a finitely generated semi-group of free group automorphisms, the log of the generalized 
spectral radius of the incidence matrices of the generators is an upper bound for the entropy 
efficiency of the semigroup. In general, finding a generalized spectral radius is difficult, but 
because of the particularly simple form of the action of the aj on the generators in Figure 2, 
in our case we can determine the exact product which yields the generalized spectral radius 
and estimate its spectral radius (Proposition 5.3). 

We conjecture that the protocol HSP^v which is used for the lower bound actually realizes 
the maximal efficiency of a point-push protocol in the presence of N obstacles when > 3. 
Figure 3(b) shows a symmetric view of the path of the single stirrer of HSP4. The path 
lies on a hypotrochoid which lead to the terminology hypotrochoid stirring protocol. The 
fact that the path of the single stirrer can be generated by a disk rolling inside a circle has 
practical implications; it implies that the stirring protocol may be implemented using a fairly 
simple system of gears to move the stirrer, see [KUlO], [FTll]. It is also interesting to note 
that the spectral radii of the matrices used in Theorem 1.1 appear to be algebraic units of 
a special type. Numerical computations for A < 20 indicate that each p{H^^'^) is a Salem 
number and each p^H^^"^) is a Pisot number. These classes of algebraic numbers arise in a 
variety of topological and dynamical circumstances, especially in optimization problems (see, 
for example, [GHOl] or [LTll]). The generator matrices whose product defines H^^^ are the 
absolute values of those which define H^^^] the connection of this fact to the algebraic type 
of their spectral radii warrants further investigation. 

Acknowledgement: We would like to thank Jean-Luc Thiffeault for suggesting this prob- 
lem and many useful conversations and thank Toby Hall for his assistance with train track 
calculations. 
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2 Preliminaries 



2.1 Linear Algebra and the generalized spectral radius 

We recall some basic facts and establish a few conventions. The identity matrix of the 
appropriate dimension in a given context is denoted /. An inequality of a pair of vectors 
means that the inequality holds in each component; the analogous convention holds for 
inequalities of matrices of the same dimensions. 

For a square matrix M, p{M) will denote its spectral radius and for 1 < p < oo, ||M||p is 
the matrix norm induced by the vector norm \\v\\p. Recall the basic facts that p{M) < \\M\\p 
and ||M||2 = (p(MM^))V2. Gelfand's formula says that p{M) = lim^^^o llM^liy"". This 
implies that if M and M' are non-negative, square matrices with M > M' , then p{M) > 
p{M'). 

For a given matrix M, we let c(M) be the row vector of its column sums, and so 
Cj{M) = Mij. An easy computation yields that for any pair of square matrices, c{MM') = 
c{M)M' . Also note that for a positive matrix M, max^(c^(M)) = 

We shall also work with free modules over the ring of all Laurent polynomials in the 
variable t with integer coefficients, Z[t^-^]. The collection of (n x n)-matrices over Z[t^^] is 
denoted Mat(n, 

We recall the definition of the generalized spectral radius. For more information see 
[Jun09]. For a finite set of matrices Ai = {Mi, . . . , M^} and a > let 

pj,[M) = max{p(M,,M,, . . . M,^f'^ : M,^ G M}. (2.1) 

Since norms are submultiplicative, 

Pk{M) < max{||M^||^ : 1 < z < n}. (2.2) 

Thus the generalized spectral radius defined as 

p{M) = lim sup /?/,(>!) 

is always finite. Note that since the spectral radius is multiplicative, for all n > 0, Pnk{M.) ^ 
Pk{M)-> and thus 

p[M) = sup{pfc(A^) : k e Z+}. (2.3) 

2.2 Word growth and the efficiency of automorphisms 

Let G be a finitely generated group with generators Xi, . . . For g ^ let £{g) be the 
minimum word length in the generators. If a : G ^ G is an automorphism, its entropy is 
defined as 

h{a) = maxlimsup ^ ^ (2.4) 

Note that h{a) does not depend on the choice of generators of G. If .4 C Aut{G) is a 
semigroup generated by {ai, . . . , a^}, for a G ^4 define its efficiency with respect to the given 
generating set as 

ett(a; |ai,...,afc|) = 



where C{a) is the minimal word length of a in the generators {a^}. Finally, the maximal 
efficiency of A with respect to the given generating set is 

Eff {ai, . . . , Qk}) = supjeff (a; {ai, . . . , ak}) : a G A}. 

Note that eff(a) and Eff(^) usually do depend on the choice of generators of A. 

Remark 2.1 If the group G is free Abelian, then an automorphism A is a linear transfor- 
mation and treating 5^ G G as a vector its word length is i{g) = ||^||i. It then follows 
that h{A) = log{p{A)). Further, as a consequence of (2.3), the efficiency of the semigroup A 
generated by a collection of linear isomorphisms {Ai, . . . , An} is the log of their generalized 
spectral radius, 

EfF(^; {A,, . . . , = log(p({A, . . . , A^})). 

Thus the entropy efficiency may be viewed as the nonlinear analog of the generalized spectral 
radius. 

2.3 Braids 

There are three of the many interpretations of braids which will be used here: as the motion of 
a collection of points in the disk, as physical braids, and as elements of a mapping class group. 
We describe each of these briefly and informally. For more information see [Bir74] , [FMar] . 
Fix a model of the two-dimensional disk and n distinguished points {pi,...,Pn} C 
Int{D'^). The disk punctured at the distinguished points is Dn := — {pi, . . . ^Pn}- 

A motion of the distinguished points is described by a family of paths rji : [0, 1] 
for i = 1, . . . ,n with r]i{t) ^ r]j{t) for all t G [0, 1] and i ^ j. To find the corresponding 
physical braid define arcs 77 : [0,1] ^ x [0,1] by fii{t) = (r/^(t),l — t). The physical 
braid corresponding to the motion of points is the image of . . . ,^n)- We can define the 
product of two point motions by following one by the other. This corresponds to stacking 
the associated physical braids one below the other. 

The elements of Artin's braid group on n-strings, denoted S^, correspond to equivalence 
classes of these objects. It can be considered as homotopy classes of point paths or else 
of physical braids. There are (n — l)-standard generators of Bn with ai corresponding to 
switching the i^^ and {i + 1)^^ strings (or points) counterclockwise. There are two types of 
relations which define the braid group: a^cr^+icr^ = cr^+icr^cr^+i for all z, and a^aj = aja^ for 

For the connection to a mapping class group, we let MCG{Dn) be the collection of isotopy 
classes of homeomorphisms of D^, where the homeomorphisms and isotopies are required to 
fix the boundary of the disk pointwise. The collection of mapping classes MCG(i^^) becomes 
a group under composition. For each generator of S^, let fi be a homeomorphism that 
switches pi and pi^i counterclockwise. The map [fi] G MCG(D^) then yields an 

isomorphism Bn = MCG(D^). 

The pure braid group on n strings is the subgroup Pn of Bn corresponding to paths rjj 
with r]j{0) = r]j{l). Thus the corresponding physical braids do not permute their endpoints 
and the mapping classes in Dn fix all the punctures. 
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Conventions: There are many possible conventions for the braid group and its various 

interpretations. Since there is not much standardization, we summarize our conventions. 
Braid words are written from left to right, and so aia2 means perform ai and then a2- In the 
corresponding physical braid, we draw (72 below ai. We maintain the same conventions when 
treating braids as mapping classes, and so aia2 means switch the first two punctures from 
the left counterclockwise and then switch the second two punctures. Thus the convention is 
the opposite of what one would use when composing functions. To be consistent with these 
conventions, when we have a matrix representation of the braid group, the matrices will act 
on a row vector from the right. Thus the matrices are the transpose of the most common 
convention. Specifically, the components of the image of the z*^-basis vector gives the i^^-vow 
of the matrix. 

2.4 Thurston-Nielsen theory and the entropy of a mapping class 

Thurston's theorem classifies mapping classes on a surface M into three classes: pseudo- 
Anosov, finite order and reducible (see [Thu88], [FLP79], [CB88], [FMar]). In each case the 
class contains what is now called the Thurston representative which has the same name, so 
pseudo-Anosov classes contain pseudo-Anosov homeomorphisms, etc. We define the (topo- 
logical) entropy of a surface mapping class x 

hix) = inf{ : / e x}- (2.5) 

Part of Thurston's classification theorem says that this infimum is achieved by the Thurston 
representative in the mapping class. In addition, the exponential word growth under iteration 
by the mapping class on 7ri(M) is equal to h{x)' Specifically, if M is a surface with chosen 
basepoint Q G M, pick a homeomorphism / G x with f{Q) = Q and let /j denote its action 
on 7ri(M, Q). Then treating as an automorphism of 7ri(M, Q) and defining h{f^) as in 
(2.4), Thurston's theorem says that /i(/j) = h{x)' For a pseudo-Anosov class the number 
A = exp{h{x)) > 1 is called the dilation 

We call a braid /? pseudo-Anosov, finite order, or reducible according to whether its 
corresponding mapping class has that type, and we define its entropy h{/3) as in (2.5), again 
using its corresponding mapping class. 

3 Point-push protocols and their topological efficiency 

In accordance with the motivation from fiuid mechanics we will sometimes use the terminol- 
ogy "stirring protocol" to refer to a braid or equivalently, a mapping class on the punctured 
disk. We study a special class of stirring protocols here in which just one stirrer moves. The 
other N stirrers (or punctures) are stationary and will be referred to as the obstacles. The 
initial and final point of the path of the single stirrer is denoted P, and so the protocol is 
given by a path F : [0, 1] ^ Dat with r(0) = r(l) = P. We may extend F to an isotopy ft on 
Dn so that ft{P) = F(t) for all t G [0, 1]. We further puncture at P forming Dat+i, and 
class of /i in MCG{Dn^i) is the point-push class corresponding to F. Since F is a closed loop 
starting and ending at P it gives an element of the fundamental group with basepoint P, 
^i{Dn^ P)' A theorem of Kra [KraSl] and Birman [Bir74] says that the point-push mapping 
class is uniquely determined by this element in the fundamental group and so elements of 
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this class of stirring motions are sometimes called TTi-protocols. 



3.1 The subgroup of point-push protocols 

Each point-push protocol with obstacles yields a mapping class that fixes the punctures 
and so they correspond to elements of the pure braid group on (A^ + 1) strings. The collection 
of such classes is a subgroup we denote Vn- If h : Pn+i Pn is the homomorphism which 
deletes the left-most string of the pure braid, then Vn = ker(/i). Equivalently, the point-push 
protocols correspond to mapping classes on D^+i with the property that when the left-most 
puncture is filled in, the resulting class in M.CG{Dn) is the identity. 



To compute the entropy efficiency we need a specified set of generators of Vn- For 
i = 1, . . . , A^, let 



where the are the standard generators of ^a^+i. Note that we have adapted the convention 
that the stirrer P is the left-most puncture and the obstacles are the next punctures. Thus 
in terms of a stirring motion, the generator ai corresponds to the stirrer passing in front of 
the first (i — l)-obstacles from the left, going around the i^^ obstacle (which is the i + 
puncture) once counterclockwise, and then returning to its initial position passing once again 
in front of the first {i — l)-obstacles (see Figure 1). The generators aj are the members of 
the standard generating set of the pure braid group Pn+i in which only the left-most strand 
is nontrivial. The next proposition follows from the theorem of Kra [KraSl] and Birman 
[Bir74] described above. 

Proposition 3.1 The group Vn of point-push protocols with N obstacles is a free group on 
the generators {o^i, . . . , Un} and is naturally isomorphic to tii{Dn-> P)- 



3.2 An adapted set of generators of tii{Dn+i) 

In working with point-push protocols with obstacles it will be frequently useful to use a 
collection of loops in 7ri(Z?Ar+i) which generate and are adapted to the protocols in the sense 




(3.1) 




Figure 1: A^" = 4: Left, the path of as; Right, the braid of 
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that the generators aj of Vn act particularly simply on them. The case = 4 is shown in 
Figure 2. The square represents the stirrer P. Fix a basepoint Q, and for i = 1, . . . , A^, the 
loop 6i begins at Q and goes counterclockwise around the i^^ obstacle. The last generator 
Sn^i encloses the stirrer and all the obstacles going clockwise around the disk. 





Figure 2: Left, the adapted generators 5j for 7ri(i^5). The arcs are used in the construction 
of the covering spaces in §4.3. Right, the action of a^. 



The reason for the last convention is to get cleaner formulas for the action of the aj on 
the 6i. We may assume that the generators {o^i, . . . , a^} of Vn act while fixing the basepoint 
Q. Figure 2 right shows the action of when = 4. Writing Q^i(5j) for the action of on 
the loop we have 

^j(^j) = (^7-1^7-2 • • • ^1^^n\i^N^ • • • • • • ^nSn^iSi . . . Sj-2Sj-i) 

ai{5j) = 5j for j ^ i. 
3.3 Topological efficiency of point-push protocols 

For a point-push protocol x ^ pick a homeomorphism h G x with h{Q) = Q, and let 
be the induced automorphism of 7ri(i^Ar+i, Q) = F^+i- Now eS{xf ^ - - - ^ ^n^}) is 
independent of the choice of h and and we use the notation 

eS{x) := efr(xa; . . . , a%^}). 

The maximal efficiency over the group Vn of point-push protocols with obstacles is denoted 
here as 

Eff (iV) := sup{eff(x) : X e Vn}- 

Remark 3.2 To avoid potential confusion we emphasize that a point-push protocol with 
obstacles is determined by an element of the fundamental group of the A^-punctured 
disk, 7Ti{Dn)^ but after also puncturing at the stirrer P it yields a mapping class on the 
(A^ + l)-punctured disk or a braid on (A^ + l)-strings. 
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4 The lower bound 

Since we are seeking Eff(A^), the maximal efficiency of a point-push protocol with ob- 
stacles, the efficiency of any single x ^ gives a lower bound. In the first subsection we 
define HSPat G Vn which numerically appears to give the maximal efficiency, and we shall 
see that its topological efficiency converges as ^ oc to the upper bound given in §5. Most 
this section is concerned with the computation of the topological entropy of HSPat- 

4.1 The hypertrophied stirring protocol 

Using the generators ai of Vn we define the hypotrochoid stirring protocol with obstacles 
as 

HSPat := Q^iQ^2 • • • o^N. 

As with any braid (see §2.3), in various contexts we treat HSFat as a physical braid, a motion 
of a collection of points (in this case just the single stirrer is moving), and a mapping class. 
Figure 3 left shows the corresponding motion of the stirrer when A^ = 4. Figure 3 right 
gives a topologically equivalent but more symmetric view in which the path of the stirrer is 
a hypotrochoid. 




Figure 3: HSP4: Left, in the braid configuration; Right, on a hypotrochoid 



The train track of a mapping class is a collapsed version the unstable foliation (see 
[PH92], [CB88], [FMar]). Given combinatorial data about a mapping class, it is always 
computable in finite time by the Bestvina-Handel algorithm ([BH95]). The edges of the 
train track give a Markov partition for a pseudo-Anosov representative in the class, and 
the induced transition on the edges from the action of gives the corresponding transition 
matrix. Figure 4 shows the invariant train track for HSP4. Its invariance under the action 
of HSP4 as well as the fact that its transition matrix is strictly positive can be confirmed 
by a rather tedious and long graphical calculation or with a computer implementation of 
the Bestvina-Handel algorithm ([Hal]). For other values of A^ the repeated regularity of 
HSPat ensures that the train tracks in those cases is constructed by a repetition of the 
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simple one-prongs connected by a single edge seen on the last three right-hand punctures of 
HSP4. The track demonstrates the important fact that the only singularities of the invariant 
foliations of HSPat are one-prongs at the punctures and a singularity on the outer boundary. 
Also note that the train track and its transition matrix prove that each HSPat is a pseudo- 
Anosov mapping class, though this is verifiable by elementary means or from Kra's Theorem 
([KraSl]). 




Figure 4: The invariant train track for HSP4. 



4.2 The method for computing /i(HSPAr) 

In principle for any one can use the invariant train track of HSPat to compute a Markov 
transition matrix for (j). The dilation of HSPat is then the spectral radius of this matrix. 
However, this computation appears intractable and so we resort to homological methods. 

As noted in §2.4, the topological entropy of a pseudo-Anosov map can be computed as 
the exponential growth rate of words in the fundamental group under the induced action of (j). 
In practise, this computation is very difficult. The corresponding computation in homology 
is much simpler, but usually just gives a lower bound ([Man74]). Fried in [Fri86] develops 
the idea that the action of the induced map on homology when lifted to finite covers (or 
equivalent ly, on twisted homology) is usually computable, often gives a better lower bound, 
and in certain cases, the bound is sharp. 

To implement this program, rather than working with individual finite covers it is simpler 
to construct a Z-covering space M and study its finite quotients. The first homology of M 
is treated as a module over with t representing the generator of the deck group. If the 

given map / (or the mapping class) lifts to M and commutes with the deck transformation 
T, then the lifted action of / is given by a matrix B{t) G Mat(m, t~-^]), where m is 
the rank of Hi{M) as a module. One may then show (see Lemma 3.2 in [BB07]) 

that a substitution of t = exp(27ri]9/g) into B{t) gives the action of / on a piece of the first 
homology of the Z^-cover constructed by taking the quotient of M by T^. 

In the case of the n-punctured disk, there is a natural Z-cover to which every mapping 
class lifts and commutes with the deck transformations. The matrix B{t) giving the action 
of the mapping class (or braid) is the Burau representation and indeed, this is now the most 
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common descriptions of the topological meaning of the Burau representation. Computing 

entropy of braids using substitutions into the Burau representation is considered in [BB07], 
and that paper is a good source for more information on the techniques used in this section. 

4.3 The Z-covering spaces 

The mapping classes representing point-push protocols are rather special and so the general 
cover used in Burau representation is not the best way to proceed. Instead we construct 
a pair of Z-covering spaces, and Zj^, which are adapted to the particular nature of the 
action of point-push protocol. Both and are covering spaces of the (A/' + l)-punctured 
disk, Dtv+i. The constructions are based on the generating loops {5j} for 7ri(DAr+i) given 
in §3.2. We shall work with simplicial one-chains and give the relation to homology in 
Proposition 4.2. 




Figure 5: Top, the spine of the covering space Z4; bottom, the spine of the covering space 
Z4. Only selected simplices are labeled and have their orientations shown. 

Geometrically, both Z-covering spaces are constructed by cutting along arcs "dual" to 
the generating loops 5i, . . . , ^at+i of 7ri(DAr+i, Q)- The arcs 7^ for j = 1, . . . , A/' connect the 
j^^ obstacle to the stirrer, and 7Ar+i connects the stirrer to the outer boundary (see Figure 2). 
For the cover Z^ one cuts along all these arcs and for one cuts along all of the arcs except 
7Ar+i. After cutting along the arcs, copies of resulting space are then stacked and pasted 
into a Z-cover. 

For algebraic descriptions, recall that a covering space with Abelian deck group of a 
manifold M corresponds to a subgroup H C Hi{M; Z), and the deck group of the resulting 
cover is Hi{M;Z)/ H . Thus a Z-covering space of D^^i is given by an epimorphism l : 
Hi[Dnj^i] Z) Z. To define this epimorphism for our two covering spaces, for j = 1, . . . , A^+ 
1, we let be the element of Hi{Dn+i]'^) corresponding to the generating loop Sj. Thus 
the set {Ai, . . . , A^v+i} gives a basis for i/i(D7v+i; Z). The covering spaces Z^ and Z^ are 
then respectively defined by epimorphisms 

i{Y^ CiAi) ^^Ci and l{Y^ qA^) = c^. (4.1) 
1=1 1=1 1=1 1=1 
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In both cases the generator of the deck group will be denoted T. 

We now give the description of the one-chains in as a free module. The wedge 

of (A^ + l)-circles given by the union of the generating loops {Sj} of 7Ti{Dn^i, Q) is denoted 
Ajv and it is a spine of D^^i. We work with the simplicial homology of the one-skeleton A^v 
and its hft to Z^, denoted Ajsf. The simplicial one-chains over Z are denoted Ci{An) and 
Ci{A]\i)^ respectively. Fix a lift Q of the basepoint Q. For j — 1, . . . , A^ + 1, let hij be the 
oriented arc (or one-chain) in Z^ which is a lift of Si and connects Q to T{Q). The collection 
{ki^ . . . , /^AT+i} will be our basis for Ci{Kn)^ and for n G Z, the formal expression f^K.j refers 
to the one-chain T^^Ki)^ where : Ci{Kn] Z) Ci{K]s[] Z) is the induced action of the deck 
transformation T. The one-chains in Aat are thus the module Ci{Zn] Z) = Z[t^^]^+^, where 
the vector of Laurent polynomials (gi, . . . , qN+i) with each qj{t) G Z[t, t~^] is identified with 
the one-chain in Ajv given by ^(lj{T^)i^j' 

The one-chains in the spine A'^ of Z'^ have a similar description. However in this case 
(^AT+i doesn't lift to an arc but rather lifts to a loop based at Q which is denoted /^tv+i- 

We 

then use the simplicial one-chains {/^i, . . . , /^a^, ^tv+iI ^ basis for Ci{K'j^) over Z[t^^] and 
proceed as with Z^- 

4A The lifted action of a point-push protocol 

We next consider the lifted action of a point-push protocol x ^ T^n- Pick a homeomorphism 
h ^ X which fixes the basepoint Q and let be its induced action on the spine A^v- More 
precisely, let : Aat ^ Aat be a simplicial map homotopic to ro/i|^^, where r : Dat+i 
is a deformation retract. Since h fixes all the punctures in Dat+i, it acts as the identity on 
Hi{Dn^i]Z). Thus h lifts to Z^ and commutes with all deck transformations. Let h be 
the lift with h{Q) = Q. The simplicial map also lifts to h"^ \ An ^ A^. Since 
commutes with the deck transformation, : Ci{An) Ci{An) can be represented by a 
matrix C{t) G Mat(A^+ 1, Zft^'^^]). The assignment x ^ C{t) is a representation of the group 
of point-push protocols with A^ obstacles which we denote $ : Vn Mat(A^ + 1, Z[t, t~^]). It 
is clear that the construction does not depend on the choice of /i G x or its retract as long 
as h^{Q) = Q. If h^{Q) = T^{Q)^ then the corresponding matrix is multiplied by t^. Thus 
we adopt the notation x* — with the understanding that it is always computed using 
an with h^{Q) — Q. The same considerations for the cover Z'^ yields a representation 
^' :Vn ^ Mat(A^ + l,Z[t^^]), and we sometimes write = C'{t) as an abbreviated 

notation. 

Using their action on the generators 5j it is straightforward to now compute the action of 
the lift of generators aj of Vn^i on the chosen basis of simplicial one-chains in the Z-covering 
space. This procedure is described in [Fri86]. We let denote the lifted action of aj on 

the simplex Hij. 

Example 4.1 We illustrate the general result by computing the lifted action of 0^3 on A4. 
The action on tti from (3.2) is 

^3(53) = S^'S^'S^'5^'5s5^5,5i52, (4.2) 

and (ys{Sj) — 5j for j ^ 3. For the lifted action to the generators of the simplicial chains 
in A, it is immediate that &3{i^j) = i^j for j 7^ 3, so we need only compute 33(^^:3). Since 
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our convention is to represent mapping classes with homeomorphisms whose hfts fix the 
image of must be an arc beginning at Q. Since from (4.2) 0^3 ((^3) begins with 82^^ the 
beginning of 0:3 (/^3) must traverse a copy of 1^2 in the negative direction and thus must in 
fact be —t~^K2' The image of Ss continues with S^^ and so the image of continues with 
— Continuing in this manner the next stage of the image 0^3 (/^s) is always forced by 
the next generator (or its inverse) in the image 0^3(53). The full image is then 

5^3(^3) = —t~^l^2 — t~^l^l — t~^l^5 — i~^/^4 + t~^Ks + t~^l^4 + t~'^l^5 + t~^l^l + 1^2 
= K.2+ t~^{l^l - 1^2) + t~^{l^b - 1^1) + i"^(^4 - ^5) + i"^(^3 - ^4) 

The process of deriving the formulas for the action of 0^3 lifted to A4 is similar, but now 
instead of traversing the arcs corresponding to lifts of ^5 to A4, one goes around the loops 
t^Kg which are the hfts of 5^ to A4. 

The same procedure in general yields the following. 

(5i(a^i) = Kn+I + t~^{i^N - i^N+i) H h t"^(/^i - ^2) 

= + t~^{l^j-2 - l^j-l) H h t"^^"^^(/^i - K.2) + t~^^~^\l^N+l - ^1)+ 

t~^{i^N — i^N+i) + t~^^^^\K.N-i — i^n) + ' • — i^j+i) when j > I 

OLi^Kj) — Kj when j ^ i 

Letting a[ the lifted action of ai on one-chains in A' we have in general 

a[[K.i) = + K.N+1 + t~^{l^N - 1^1 - I^N+l) + t~'^{l^N-l " I^n) H h t~^^~^\hi2 - l^s) 

&2{'^2) = i^AT+l + {l^N - ^AT+l) + t~^{l^N-l " I^n) H h t~^^~^\Ki - K2) 

= + t~^[Kj-2 - f^j-l) H \- t~^^~^\K2 - t^z) + t~^^~'^\Ki - K2 + + 

— /^i — K.N^i) + t~^{f^N-i — ^at) H h t~^^~^\Kj — /^j+i) when j > 1 

— K.j when j i 

(4.4) 

4.5 The double cover 

The two-fold covering spaces formed by taking the quotient of and by the square 
of the deck transformation will be central to our calculations. We recall some basic facts 
concerning the connection of various objects in the full Z-covering spaces, the base space 
and this two- fold cover, and then compute the action on homology in this cover. We will 
initially just discuss and comment on afterwards. To simplify notation fix N and 
drop it from the notation, and so now Z = Z^. 

Recall that the Z-covering space p : Z ^ ^at+i corresponds to an epimorphism l : 
Hi{Dn^i;Z) Z. The two-fold covering space p^^^ : Z^^^ ^at+i corresponds to the 
epimorphism z.^^) ; Hi{Dn+i] Z) Z2 where z.^^) ;= 7^(2) o with tt^^) : Z ^ Z2 = Z/(2Z) the 
projection. Thus Z^^^ = Z /T'^, and we denote the deck transformation of Z^^^ as T^^^. The 
corresponding one-skeleton of Z^'^^ is denoted A^^^. 



14 



If T is an indeterminate with = 1, then Ci(A(^)) is a Z[r] module with Ci(A(^)) = 
Z[r]^"^^. The action of ai is denoted as af^ and using (4.3) we get 

(-1^(1 - r)^i + {-iy^\l - t)k2 + ••• + (!- r)/^,_i + k,- 

(1 - t)k^^y H V - r)/^7v+i when is even 

(-l)^'(l - r)/^i + (-iy+'(l - t)k2 + ••• + (!- T)K,_r - TKj+ (4.5) 
(1 - t)i^j+i + • • • + (-l)(^'+^Hl - r)i^N+i when is odd 
tij for all N when j 7^ i 

Similarly, Z'^^) = Z' /T'^ is defined by the epimorphism l^'^^ :— tt^^^ o l\ Its one-skeleton is 
A'(2) and Ci(A'(2)) ^ Z[r]^+i and using (4.4) the action of 5^ is 

• + (1 - r)K.j_i - TKj + 

when A^ is even 

' + {l-T)tij_i + i^j- (4.6) 
when A^ is odd 

Let X be a point-push protocol treated as a mapping class. The projection of the action 
X* on Ci(A) to one on Ci(A(2)) is denoted xf^. The action of x on i/i(A(2); Z) ^ i/i(Z(2); Z) 
is written The same notations are used for Z'^'^\ but with a prime attached. 

The next proposition gives the connection of the substitution t = — 1 into the matrices 
C{t) and C\t) and the growth rate of the action of x on the first homology of the two- fold 
cover Z^'^\ It follows from Lemma 3.2 in [BB07] and the fact that the growth on chains 
is the same as that on homology. For sake of having the paper self-contained, we give the 
elementary proof. 

Proposition 4.2 Assume that x ^ T^n is a given point-push protocol As constructed above, 
Z^^) and Z'^'^^ are double covers and x^^ and Xi^'^^ the induced actions on i/i(Z(^^) and 
Hi^Z'^'^^), respectively, and $ and are representations o/Vn- Letting C{t) — $(x) ci'^d 
C\t) = then 

Pix^^^) = PiCi-1)) and p(xf )) = piC'i-1)). 

Proof: We will initially just discuss Z^ and comment on afterwards. 

Since (T^^^)^ = id, the eigenvalues of its action on simphcial one-chains in A^^^ are ±1, 
and we denote the corresponding eigen-space decomposition as Ci(A^^)) = ® The 
subspaces V± are also characterized as V+ = {c+T^^\c) : c e Ci(A(2))} and V- = {c-T^ : 
c e Ci(A(2))}. Since xl^^ commutes with tP , we have that x*^ preserves the decomposition 

Ci(A(2)) = wey- 

Treating Ci(A'^^)) = Z[t]-^"'"^, the eigen-subspaces have alternative descriptions as = 
{c + TC : c E Ci(A(^))} and V- = {c — re : c E Ci(A(^))}. The homomorphism x^*^ acts on 
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d'PiK,) = (-1^(1 - r)Ki + (-ly+Hl - t)k2 + ■ ■ 

(1 - T)Kj+i + ■■■ + {-iy{l - t)kn+1 

= (-1^(1 - + {-iy+\l - t)k2 + ■ ■ 
(1 - r)Kj+i + ■■■ + (-1)^(1 - t)kn+i 
af^\K.j) = K.j for all A'^ when j i 



Ci(A(^)) as the matrix C(r), and we will write C{t) = Bq + tBi with Bi G Mat(7V + 1,Z). 
Thus acting on an element of we have (c + rc)C(r) = {cBq + cBi) + t{cBi + cBq) — 
{c+Tc){Bo + Bi) = {c+Tc)C{l) Similarly acting on an element of V- we have {c — tc)C{t) = 
(c — Tc)C{—l). So xl^^ acts on V+ by C{1) and on V- by 1). Substituting r = 1 into 
(4.5), we find C(l) = / (this also follows from the fact that pl^^ : Ci(A(2)) ^ Ci(D7v+i) 
is an isomorphism when restricted to V^). Since det(C(— 1)) = 1, p{C{—l)) > 1 and so 

To find out what this means for homology we go back to working over Z. For 1 < i < 
A^ + 1, let l^i = Ti'^\tvi)^ and so {Hii^l^i} forms a basis for (7i(A(^);Z). Now for 1 < i < A, 
let Ti = tvi — tvi^i and F^ = — l^i^i. Letting F = k^n+i + ^aa+i it is easy to check that 
{Fi, . . . , Ftv, Fi, . . . , Fjv, F, i^N+i} is also a basis for (7i(A(^^; Z) and we define V as the linear 
spanj/ = (Fi, . . . , Fat, Fi, . . . , F^v, f ) Thus Ci(A(^); Z) = V®{kn+i)' Now f and all of the F^ 
and F^ are cycles because there are no two cells in A^^^ In fact, V is a concrete realization of 
the simplicial homology i7i(A(^^; Z). Either by direct calculation using (4.5), or from general 
principles, x^\v) = V. Also from (4.5), a^^\KN^i) = /^at+i, and so we see that the spectral 
radius of xl^^ acting on all of (7i(A(^^; Z) is the same as the spectral radius of its restriction 
to V. But xl^^ restricted to V is exactly the induced action of /i^^^ on i7i(Z(^);Z) and 
so by the conclusion of the previous paragraph, p(xi^^) = C{—1). 

For the other cover the proof on the chain level that p(x*^^) = P{C\~^)) is exactly 
the same as for Z^'^\ A slight alteration is required when passing to homology because ^at+i 
lifts to a loop in Z'^^^: for 1 < i < A — 1, let F^ = /^^ — and F^ = "^^ — "^^+1, and 
finally let T = k.i It is easy to check that {Fi, . . . , Fat-i, i^n+i^ Ti, • • • , Tat, ^a^+i, ^i} 

is a basis for (7i(A'(^); Z), and we let = (Fi, . . . , Fat-i, /^a^+i, Fi, . . . , Fat, /^a^+i, F). Thus 
Ci{A'^^^]Z) = V'® (/^i), and V is a concrete realization of Hi(k'^^^',Z) with x^'^\v') = V . 

f(2) /(2) 

From (4.6), (i^i) = r/^i + . . . when A^ is even, and (i^i) = /^i + . . . when A^ is odd, 
and aj^\Hii) — Hii for j 7^ 1, and so the spectral radius of x*^^ acting on all of Ci(A'(^);Z) 

is the same as the spectral radius of its restriction to V\ But recall that x*^^ restricted to 
V is exactly x^'^\ the induced action of x'^^^ on i^i ( A'^^^ ; Z) . Therefore by the conclusion of 
the first paragraph, pix^'^^) = C''(— 1). ■ 

4.6 Another representation 

Proposition 4.2 shows the significance of the substitution t = — 1 into the matrices C{t) 
and C\t). It will turn out that the matrix of significance for computing /i(HSPAr) will be 
C{t) when A^ is even and C'{t) when A^ is odd. Conveniently, the substitution of t = — 1 
into (4.3) (or r = — 1 into (4.5)) when A is even yields an identical formula to the result of 
substituting t = — 1 into (4.4) when A is odd. 

Now note that for all j, aJ'\K.N+i) = /^at+i, and so for any product of the af \ at is an 
eigenvector with eigenvalue 1. Since we are interested in exponential growth rates, we may 
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ignore k^n+i- For 1 < k < N, define the (A^ x A^)-matrix T^^^ by 

T^f = 2 . (-1)^-^+1 for 1 < j < A; - 1 

T^l^ = (4.7) 

Tj^f = 2 . (-1)^-^ for A; + 1 < j < 

T^-^ =0 for z ^ fc, 



and let 



Extend the assignment ak ^ E^^^ to a homomorphism : Vn SL(A^, Z). Thus ^^(x) is 
the matrix C{—1) with the last row and column deleted for N even and ^^(x) is the matrix 
C(— 1)) with the last row and column deleted for N odd. 

4.7 A trace computation 

Fix an > 2. Recall that HSPat — ai . . . a^. We will compute trace(^^(HSPAA)) inductively 
using the initial segments of the product ai . . .ak for k = 1, . . . , A^. These products have 
the matrix representation H^^^ := . . . o^fc) = E^^^E^'^^ . . . E^^\ 

Lemma 4.3 With * and HSPtv as defined above, trace(^^(HSP7v)) = -3^ + 3A^ + 1. 

Proof: It follows easily from the definitions that for any 1 < a, 6 < A^, 

We now prove by induction on k < N that for m = 1, . . . , A^ — /c, 

^SVm = 2 • 3'=-X-l)^+'=+"^ 1 < i < A; 

i^^iVm = A; + l<i<A^andi7^/c + m (4.9) 

i/^'^^ ' = 1 

The case A; = 1 is trivial. So assume (4.9) is true for k. Thus using (4.8). for 1 < z < A; and 
m = 2, . . . , A^ - A:, 

tr(/c+l) _ rrik) , trW 7.(^+1) 
^i,k-\-m ^i,k+m /c+l,/c+m 

2 . 3^~^( 2)^"'"^"'""^ -|- 2 • 3^~^( '^Y~^^~^^ • 2( 2)^~'~"^~(^~'~''"'^ 

as required. The other values of i are special cases easily checked, completing the proof of 
(4.9). 
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Now we prove that for 1 < A; < iV, trace(i/('^)) = -3'' + 2k + N + 1, again by induction 
on fc, with the case k — 1 being trivial. Using (4.8) and the formulas (4.9), 

N 

trace(i7(^+^)) = trace(i/(^)) + J2<LTiltl 

a=l 
k 

= trace(FW) + • 3fc-«(_i)«+fc+i2(-l)^-« 

a=l 

^_3^+i + 2(A; + l) + A^ + l. 
The statement in the lemma now follows by letting k — N . ■ 

4.8 The orientation double cover of a pseudo-Anosov map 

A pseudo-Anosov map (j) is characterized by the existence of a pair of transverse invariant 

measured foliations. If these foliations are oriented, it is now standard that the dilation 

of the pseudo-Anosov map is the spectral radius of the action of (j) on first homology. For 

pseudo-Anosov maps on the punctured disk and in many other cases, a pseudo-Anosov map's 

foliations are not oriented. However, one can lift (j) and its foliations to a new pseudo-Anosov 
— * — * 

map on a branched cover M where there foliations are oriented, and thus the dilation is 
detectable using homology. The construction of this cover requires a fair amount of a priori 
information about 0, and the exact connection between (f) and (j) can be difficult to ascertain, 
so in many situations the orientation double cover is not a useful tool. However, because the 
family HSPat lives on a genus zero surface and its foliations are of a particularly simple type, 
it will turn out that its orientation double cover always corresponds to one of the two-fold 
covers finite covers constructed by taking quotients of the Z-covering spaces Z^'^^ and Z'^'^^ 
which were considered in §4.5. 

A few definitions are required to give a characterization applicable to the case at hand. 
Assume is a pseudo-Anosov map on the n-punctured disk with unstable foliation T'^. 
Let S be the collection of singularities of T'^ which are of odd order and contained in the 
interior of D^. In addition, P denotes the collection of punctures of and S' d P the 
subset consisting of those punctures at which the singularity of J^^ is of odd order. The 
orientation double cover of (j) is the two- fold cover branched over S \J S' . More formally, 
first further puncture at S to form M' — — S. For each pi G P U let F^ be 
a small loop going counterclockwise around p^, and so the homology classes {[F^]} are a 
basis for Hi{M'] Z). If / is the set of indices i so that pi ^ S' VJ 5, define the epimorphism 
Hi{M']Z)^Z2 by 

Q[Fi] ^ ^ Q mod 2. 

The corresponding two- fold cover over M' is denoted M^ Finally, fill in the punctures in M' 
that project to points in S to form M" which is a branched covering space over D^. The 
following is standard: for example, it follows from Lemma 2.2 in [BB07]. 

Proposition 4.4 Let (j) he pseudo-Anosov map on the n-punctured disk and M" the 
branched double cover of formed by branching over the singularities of the foliations of 
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(/) which have odd order as was just constructed. Then M" is the orientation double cover of 

0. 



4.9 The entropy of HSPat and the lower bound 

We now have all the ingredients needed to obtain our lower bound. 

Theorem 4.5 The topological entropy of the hypotrochoid stirring protocol with N obstacles 
is given by the log of the spectral radius of the matrix defined in the previous subsection, 

/.(HSPat) = log(p(i/(^))). (4.10) 

As a consequence, 

log( ) < HRSPm), 

and thus the maximal efficiency of a point-push protocol with N obstacles satisfies 

log(3^ - 3A^ - 1) - log(A^) 



< Eff(7V). 



Proof: The main observation we need is that the double covers Z^^^ and Z'^^^ from §4.5 
are identical to the orientation double covers of a pseudo-Anosov map G HSPat when 
N is even and odd, respectively. Once again we discuss the case of Z^^^ first. Recall that 
2^(2) corresponds to the epimorphism p^^^ : Hi{Dn^i]Z) Z2 given by P^'^\Y^fJi^ CiAi) = 
^i=i^ Q mod 2 with the basis given in §4.3. Let Q be the homology class of a small 
loop going once counterclockwise around the stirrer P. We then have Q = — Y2i=i^ 
so p^'^\^) = 1 mod 2 when is even. As noted in §4.1, the foliations of the pseudo- 
Anosov representative (j) G HSPtv have their only singularities at the punctures and outside 
boundary, and the singularities at the punctures all have odd order. Thus by Proposition 4.4, 
the orientation double cover of (j) is defined by the epimorphism, p(A^) = 1 mod 2 for 
i = 1, . . . , A^ and p{^) = 1 mod 2. Thus the orientation double cover of HSPat is exactly 
when A^ IS even. 

Again for even A^, letting C{t) = ^{HSPn), Proposition 4.2 says thatp(0i'))=p(C(-l)), 
where 0^^) is the hft of the pseudo-Anosov(/) G HSPat to Z^'^\ But since Z^'^^ is also the orien- 
tation double cover for 0, we have h{(f)^'^^) = log(p((^i^^)). Since entropy is preserved by finite 
covers, we also have = log(p(0P)). Finally, as remarked in §4.6, p{H^^^) = 
finishing the proof when A^ is even. 

To obtain the result when N is odd, we use the two- fold cover Z'^'^K It corresponds to the 
epimorphism /^)(Ez=t^ ^i^i) = E^Ii mod 2. Thus when A^ is odd, p<^\n) = 1 mod 2, 
and once again Z'^^^ is the orientation double cover of the pseudo-Anosov G HSP^v. The 
rest of the argument yielding (4.10) is the same as when A^ is even. 

Since for any (A^ x A^)-matrix C, p{C) > trace(C)/A^, the last two statements of the 
theorem follow from Lemma 4.3. ■ 
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4.10 The case = 2 

The case of two obstacles requires separate consideration. While Theorem 4.5 is true for 
N = 2^ the lower bound is simply < /i(HSP2) and indeed, it is easy to check that 
p(HSP2) = 1. Fortunately, there are well-known, simple methods for analyzing the required 
mapping classes. The Euler-Poincare formula ([FLP79]) implies that any pseudo-Anosov on 
the thrice-punctured disk has invariant foliations whose only singularities are one-prongs at 
the punctures and a one-prong on the outer boundary. As in the proof of Theorem 4.5, this 
implies that for any point-push protocol x ^ ^2, h{x) — log(p(^(x)) with ^(x) defined in 

§4.6. 

For j = 1,2, let Nj = "^{o^j). Now x is a word in the generators and thus 

with exactly the same indices, ^(x) is a product of matrices from the set Af :— {A^i^^, ^2"^}- 
Finding the maximal efficiency in V2 is then exactly the same as finding the generalized 
spectral radius of Af (see Remark 2.1). From (4.7), 

2-(-l)±^\ 1 0\ 

-[0 1 J' - 1^2- (-1)^1 ij- 

The following proposition is no doubt well known; we include the simple proof. 

Proposition 4.6 If A4 = {Mi, . . . , M^} and the matrix Mi has the property that ||Mi||2 = 
max{||M^||2 : 1 < i < n} and further^ G Ai^ then the generalized spectral radius of the 
collection M is p{M) = p2{M) = \\Mi\\2 = {p{MiM^)f''^ . 

Proof: By (2.2), p{M) < \\M\\2, and since pi^MiM^f'^ = l|M||2, the upper bound is 
realized by MiMf . ■ 

Theorem 4.7 The maximal entropy efficiency for point-push protocol with two obstacles is 
realized by the protocol olxOl^^ ^ and so has value 

Proof: The set M is closed under taking transposes and all its elements have the same 
two-norm, namely (3 + 2v^)^/^ = 1 + a/2, and so the theorem follows from Proposition 4.6. 
■ 

The same argument obviously also shows that the maximal efficiency is also realized by 
a{^a2 or a20L{^ or a^^ax. Also notice that in terms of the usual generators of the braid group 
on three strings, a\a2^ = g\g\G2^Gy^ . Conjugating by g\ we get g\g2^^ and so the path 
of the stirrer is a figure eight with an obstacle inside each loop. This protocol is called the 
"taffy puller" because it is used in making certain kinds of candy (c/. Section II in [FTll]). 

5 The upper bound 

The main tool for finding an upper bound to the entropy efficiency is the incidence matrix of 
the action of a point-push protocol on the fundamental group of the punctured disk. Propo- 
sition 5.1 says that the generalized spectral radius of the collection of incidence matrices gives 
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an upper bound for the entropy efficiency. The computation of the appropriate generahzed 
spectral radius is possible here because of the particularly simple form of the action of the 
generating automorphisms on the generating loops 5i of the free group 7ri(DAr+i, Q) given 
in (3.2). 

5.1 Incidence matrices and the generalized spectral radius 

We start with a general result that says that the generalized spectral radius of the collection 
of incidence matrices is an upper bound for the efficiency of a collection of automorphisms 
of a free group. 

Assume now that the group G — F^, the free group on the generators Xi, . . . Let 
ni(g) be the number of occurrences of the generator Xi or its inverse x^^ in the reduced 
form of g e G. We denote the row vector n{g) = {ni{g), . . . , nm{g)), and so ||n(5')||i = £{g), 
the word length of g. The incidence matrix of the automorphism a : is the 

(m X m)-matrix A with Aij — nj(a(xi)). 

Proposition 5.1 //{ai, . . . , ak} is a collection of generators of the semigroup A C Aut(F^) 
whose incidence matrices are {Ai, . . . , Aj^}, then the entropy efficiency of A is bounded above 
by the generalized spectral radius of the collection of incidence matrices or 

ES{A; {ai,...,a,}) < . . . , ^,}). 

Proof: Let a e A have incidence matrix A. By construction, n{a{g)) < n{g)A, and iterating 
we have for n > 0, 

n(a"(^)) < n{g)A^. 
Taking one-norms, n*'^ roots, and then using Gelfand's formula 

limsup ||n(a"(^))||i/" < limsup ||n(^)||i/"P"||i/" = p{A). 

Since this is true for all g G F^, using the definition in (2.4) and taking the logarithm we have 
h{a) < \og{p{A)). Now treat a as a product of the generators in reduced form, a — ai-^ . . .ai^. 
Since incidence matrices ignore cancellations, A > Ai^ . . . ^4^^, and thus because the matrices 
are positive, p{A) > p{Ai^ . . . A^^). Therefore for all a of word length /c, 

and the result follows. ■ 

5.2 The incidence matrices and the form of the upper bound 

The action of ak on the generating loops 5i of the free group 7ri(DAr+i, Q) as given in (3.2) 
has {N + 1) X {N + l)-dimensional incidence matrix A^^'^ — I + 5'^^^ with 5'^^^ defined by 
_ ^g) _ 2 for j ^ fc, and S^-^ = otherwise. As in §4.6, for all fc, the last row oiA^^^ 
is (0 ... 1). Thus if we define the reduced incidence matrix A^^) as the matrix obtained 
by deleting the last row and column of A , then 

p{{A^'\...,A(^^})=pi{A''\...,A'''^}). (5.1) 
21 



Thus we may restrict attention to the A^^\ We still have A^^^ = / + S^^\ but now S^^^ 
is N X N. We could have alternatively defined A^^^ as the absolute value of the matrix 
'^^{ak) = E^^"^ giving the action of ak on lifted chains to the two-fold cover from §4.6. 

It is easy to check that A is also the incidence matrix for the functional inverse, of 
ak and further that H^l^^^lli = 3 for all k. Thus by Proposition 5.1, (2.2) and (5.1) we have 

Proposition 5.2 If {A^^\ . . . , A^^^ is the set of reduced incidence matrices of the standard 
set of generators {af^^ . . . , a^^} of the group of point-push protocols Vn, then the maximal 
entropy efficiency in the presence of N obstacles satisfies 

Eff(7V)<p({AW,...,A(^)})<log(3), 

where p{{A^^\ . . . , A^^^) is the generalized spectral radius of the set of incidence matrices. 

5.3 Computing the generalized spectral radius 

We can improve the upper bound of log (3) in Proposition 5.2 by computing the generalized 
spectral radius p{{A^^\ . . . , A^^^})^ at least in the sense of showing that it is realized by an 
explicit matrix. 

Proposition 5.3 If H^^^ := . ..A(^\ then p{H^^^y/^ = p{{A^^\ • • -.A^^^), 

We shall see in Remark 5.7 that the product realizing the generalized spectral radius in 
Proposition 5.3 is not unique. 

The computation of the generalized spectral radius makes frequent use of the fact that 
multiplication by an incidence matrix A^^"^ changes the column sum in a fairly simple fashion. 
Recall that c(M) denotes the column sums of M. We then have 

Cj{MA^^^) = Cj{M) + 2ck{M) for j ^ k, and Ck{MA^^^) = Cfc(M). (5.2) 

Now the maximal column sum of a positive matrix is its one-norm, ||M||i = max c(M). Thus 
if we can determine a sequence of multiplications by incidence matrices which maximizes the 
column sums, then we know the sequence which maximizes the one-norm. An argument 
using Gelfand's formula then completes the proof. 

Finding the proper sequence is not difficult, again because of the simple action of the 
j[W on column sums. An example will make it clear. 

Example 5.4 We begin with the identity matrix / and c(/) = (1,1,..., 1). Since all the 
entries are equal, it doesn't matter which A^^^ we choose to multiply by, and so for simplicity 
we choose A^^\ and we get c(M(i)) = (1, 3, 3, 3, . . . , 3). Now to maximize the column sums 
in the next product, (5.2) indicates we should choose an index k for which the corresponding 
entry in c{IA^^'^) is maximal. Again, we may choose any > 1, but for simplicity we choose 
yl^^) yielding c{IA^^'^A^'^^) = (4, 3, 6, 6, ... 6). Continuing we see that at the /^-step, the 
entry of the current column sum will be among the maximal values and so we choose to 
multiply by A^^\ It is more or less clear then that after steps we will get a maximal 
column sum by using the sequence prescribed in H^^"^ := A^^^ . . . A^^^ 



22 



We need a few definitions and preliminary lemmas to formalize the insight of the example. 
Fix N > 3 and let O^v be all A/'-tuples of positive integers given in non-increasing order, so 

ON = {de (Z+)^ : ai > a2 > • • • > a^v > 0}. 

The sorting function is the map S : (Z+)^ On which sends a row vector of positive 
integers to the row vector of its entries sorted in non-increasing order. In keeping with our 
conventions we will have 5' act from the right, and so 65' is the sorted version of 6 G (Z+)^. 

It is easy to check that for each j and each a G (Z"^)^, aAj G (Z"^)^, and so we may 
define Wj : On ^ On hj aWj := aAjS. Using the fact that all entries of a are strictly 
positive we get the formulas 

aWi = (a2 + 2ai, . . . , + 2ai, ai) 

aWj — (ai + 2a . . . , aj-i + 2a a^+i + 2a an + 2a aj) when j > 1 

The next lemma collects some useful facts about the action of the Wj and the usual partial 
order on (Z+)^. Its proof is an elementary exercise using (5.3). 

Lemma 5.5 Assume a^b G (Z+)^. 

(a) Ifa>h, then aWj > bWj for all j. 
(h) If a>b, then aWj > bWj for all j. 

(c) If j < k and aj = a^, then aWj — bWk- 

(d) If j < k and aj > a^, then aWj > bWk- 

(e) If j < k and a > b, then aWj > bWk- 

The next step is to formalize the various sequences of multiplying the matrices in our 
family. A strategy is a list 5 = (si, ^2, . . . , s^) with 1 < Si < N ior all i. Each strategy s 
defines a map Ws : On On given by aWs_ := aWs^Ws^ • • • Wsj^. 

The next definition gives a canonical form for a strategy which will facilitate comparing 
the results of two strategies. Given a ^ On and a strategy 5, the standardization of s with 
respect to a is the strategy = (s^^, . . . , 5^^) defined inductively as follows. Let s[ = 
min{i : a^ = a^-^}. For m > 1, let b^^"^ = a W^/^ . . . Ws'^_^, and let 5^ = mm{i : fo-"^^ b^J2^}. 
Note that it follows easily that for all m, 6^"^) = aWg^ . . . Ws^_^^ and so in particular, 

aWs^aWs'. (5.4) 

If the strategy s is equal to its standardization with respect to a, then s is said to be standard 
for a. In light of (5.4) we henceforth restrict consideration to standard strategies. 

The simplest strategy is to repeatedly use the first place, which means using the largest 
entry; for each /c, let 1^^^ be the /c-tuple of all ones, 1^^^ = (1, 1, . . . , 1). The next proposition 
collects some useful facts about strategies and products of the matrices from M. The first 
part says that l^^^ it is always the best strategy for maximizing the end result. We adapt 
the simplified notation that given an fc-tuple £— {£i . . .£k)^ ^i — A^^^^ . . . A^^'^\ 

Lemma 5.6 Assume a e On- 
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(a) // 5 = (^i, . . . , is standard strategy for a and s ^ then aW-^ik) > aWg- 

(b) If s = {si^ . . . ^ Sk) is any strategy^ then aW-^ik) > aWg- 

(c) Given any k-tuple (-^i, . . . , -^/c); there exits a strategy s — (si, . . . , 5^) with aA^S — 
aWs. 

(d) If ^ is the k-tuple (l,2,...fc) anda^ G (Z+)^, then a A^S ^ aW^^k) . 

Proof: Since by hypothesis, s ^ there is a smallest index i with Si > 1. If i = 1, let 

h — a and if z > 1, let h = aW(^si,...,si-i) = W^{i-i){a). Since s is standard for a, when i = 1, 

as- < ai and so hg. < 6i, and when z > 1, again since s is standard for a, bg. < 65.-1 < 61. 
^ — * — * 

Thus in either case by Lemma 5.5(d), bWi > bWs- which is to say that aW-^ii) > aWs^^.^^si- 
Since Sn> 1 for all n > z, by Lemma 5.5(e) aW^(k) > aWg, as required for (a). 

For (b), if 1^^^ is the standardization of 5, then aW^ik) — aWs^^ otherwise use (a). We 
prove (c) by induction on /c, with k = 1 being trivial. The inductive step is based on the 
easily checked observation that for any (perhaps not ordered vector) d G (Z+)^, and for any 
J, if we let / be the index with {dS)jf = dj^ then dAj S = dSAjf S. 

Now assume that the result is true for a fc-tuple i. Given — (£1, . . . ,£/c,£/c+i), using 
the observation of the previous paragraph with d — aA^^ there is a / with aA^A^^^^ S — 
aAiS Ajf S = aWs_Aj' S, using the inductive hypothesis. Thus aA^^ S = aWs+ with s'^ = 
(51, . . . , Sk^j')^ completing the induction. 

The proof of (d) is a simple calculation. ■ 

Proof of Proposition 5.3: Given a A; > 0, write k = mN + j with < j < and let 

any product of k matrices 

P:= A^^^\..A^^^\ \\B^''^\i > \\P\\i. 

Let a — (1, 1, . . . , 1) e {Z'^)^. By Lemma 5.6(c), there is a strategy s which we may 
assume to be standard so that aP S — aWg. Further, by Lemma 5.6(d), a B^'^^S is produced 
by the strategy l^'^), so aB^'''>S = aW^(k). Thus by Lemma 5.6(b), a B'^''^ S >dPS. But for 
any positive matrix M, max(aM) = ||-^||i and so > ||P||i, finishing the proof of 

the first step. 

Since = (^W)m^(i) _ _ ^j^g g^-g^ g^gp implies that ||P||i < c||(^W)"^||i, where 

c = max{||^(^) . ..A^^^Wi : 1 < j < N}. 

But for any matrix M, p{M) < ||M ||i, and so for each k — mN + j letting A4 — 
we have 

PkiM) < \\{H'^''Y\\i'c'/'^ < ||(^W)'"||i/("^^)cV^ 

and so using Gelfand's formula, p{A4) < p{H^^^Y^^ . On the other hand, for any m > 
p((^W)-) = p(^W)-, and so PmN{M) > p{H^^^y/^ and thus p{M) > p{H^^^y/^ , 
finishing the proof. ■ 

Remark 5.7 The product of matrices in Proposition 5.3 is by no means unique. If i is any 
A/'-tuple which is a cyclic permutation of (1, 2, ... A/') and a — (1,1,...,!) G (Z"^)^, then 
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it is easy to check that aAiS — aW^ik). Thus the joint spectral radius of the collection 
. . . , A^^)} it is also realized by any product of N of the A^^) in which each A^^) occurs 
exactly once. 

5.4 Properties of H^^"^ and the upper bound 

Recall that Cj{M) denotes the sum of the entries in the j*^ column of the matrix M. 

Proposition 5.8 The product of incidence matrices H^^^ :— A^^^ . . . A^^^ satisfies Cj^H^^"^) — 
3^ - 2 • 3^-^ and trace(^(^)) = 3^ - - 1, and so 

_ AJ _ ^ 

^--^ < piH^""^) < 3^ - 2. 

Proof: Fix and for 1 < A; < A^, let H^''^ = A^^M^^) . . . A^'^). An easy induction yields 
that Cj{H^''^) = 3^= - 2 • 3^-^ for j <k and Cj{H^''^) = 3'' for j > k. The column sum formula 
then follows by letting k — N. 

For any matrix M, it follows from the definition of A^^"^ that trace(M74(-^^) = trace(M) + 
2cj{M) — 2Mjj. Since for any k, it is easy to check that H^^^^ = 1, and we have just seen 

that Ck{H^^~^^) = 3^~^, so another simple induction gives trace(i/(^)) = + 3^ - 2fc - 1, 
and we again let k = N . 

The estimates on the spectral radius then follow from the general fact that for any 
(A^ X A^)-matrix M, | trace(M)|/A^ < p(M) < ||M||i. ■ 

As a consequence of Propositions 5.2, 5.3, and 5.8 we have our upper bound for the 
entropy efficiency. 

Theorem 5.9 The maximal efficiency Eff (A^) of a point-push protocol in the presence of N 
obstacles satisfies 

^^^^^ N N • 
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